Prom regular to growing small-world networks 
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We propose a growing model which interpolates between one-dimensional regular lattice and small- 
world networks. The model undergoes an interesting phase transition from large to small world. 
We investigate the structural properties by both theoretical predictions and numerical simulations. 
Our growing model is a complementarity for the famous static WS network model. 
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I. INTRODUCTION 

Many real-life systems display both a high degree of 
local clustering and the small- world effect [3, H, H, 0, @ • 
Local clustering characterizes the tendency of groups of 
nodes to be all connected to each other, while the small- 
world effect describes the property that any two nodes 
in the system can be connected by relatively short paths. 
Networks with these two characteristics are called small- 
world networks. 

In the past few years, a number of models have been 
proposed to describe real-life small-world networks. The 
first and the most widely-studied model is the simple and 
attractive small- world network model of Watts and Stro- 
gatz (WS model) @, which triggered a sharp interest 
in the studies of the different pro per ties of small-world 
networks and WS model @, 0, U,M, HO, E O- The 
WS model is probably a reasonable illustration of how 
a small- world network is shaped. However, the small- 
world effect is much more general, researchers begin to 
explore other mechanisms producing small-world net- 
works. Recently, Ozik, Hunt and Ott have introduced a 
simple evolution model (OHO model) of growing small- 
world networks with geographical attachment preference, 
where new nodes are linked to geographically nearby 
ones A deterministic version of a special case of 

the OHO model was presented in and further ex- 
panded in (l6j . In addition, many authors found that 
small- world properties can be also created in other inter- 
esting ways [171, [l8| . 

It is well-known that the WS model is the first success- 
ful model that interpolates between a regular ring lattice 
and a completely random network, and plays an impor- 
tant role in network science. However, the WS model 
is static (i.e. the network size is fixed), this does not 
agree with the growth property of many real-life sys- 
tems [19]. In addition, in real systems, a series of mi- 
croscopic events shape the network evolution, including 
addition or removal of a node and addition or removal of 
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an edge [201 |21l |22l I23L |24| . Therefore, it is interesting to 
establish a growing small-world model to investigate the 
effect of local events on the topological features like the 
WS model. To our best knowledge, all previous models 
of small-world networks either are static or only took into 
account the addition of nodes and edges without consid- 
ering other microscopic events. 

In this paper, we present a growing small- world net- 
work model controlled by a tunable parameter q, where 
existing edges can be removed. By tuning parameter 
q, the model undergoes a phase transition from large to 
small worlds as the WS model. We study analytically and 
numerically the structural characteristics, all of which de- 
pend on parameter q 



II. THE MODEL 

In this section, we introduce a growing model which de- 
scribes networks from regular to small-world. The model 
is constructed in the following way (see Fig. [TJ. 

(i) Initial condition: We start from an initial state (t 

1 : of three nodes distributed on a ring, all of which form 
a triangle. 

(ii) Growth: At each increment of time, a new node 
is added which is placed in a randomly chosen internode 
interval along the ring. Then we perform the following 
two operations. 

(iii) Addition of edges: The new node is connected its 
two nearest nodes (one on either side) previously existing. 
Nearest, in this case, refers to the number of intervals 
along the ring. 

(iv) Removal of an edge: With probability q, we re- 
move the edge linking the two nearest neighbors of the 
new node. 

The growing processes are repeated until the network 
reaches the desired size. 

When q = 1, the network is reduced to the one- 
dimension ring lattice. For q — 0, no edge are deleted, 
the model coincides with a special case of the OHO 
model [14]. Varying q in the interval (0,1) allows one 
to study the crossover between the one-dimension regu- 
lar lattice and the small- world networks. 

By construction, at every step, the number of nodes 
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the form 



FIG. 1: Scheme of the growing network for the case of q — 0.3, 
showing the first four time steps of the evolving process. The 
dashed line represents the edge removed. 



increases one, while the average number of edges added 
is 2 — q. Then we can see easily at time t, the network 
consists of t + 1 nodes and average (2 — q)t + 2q — 1 edges. 
Thus when t is large, the average node degree at time t 
is equal approximately to a constant value 4 — 2q, which 
shows our network is sparse like many real-life networks. 



p(k, s, t + 1) 



2(1 ~q) 
t+1 



p{k - l,s,t) 



with the initial condition, p(k,s = 0, 1,2, t = 2) = 5k,% 
and the boundary one p(k, t, t) = 8k,2- This accounts for 
two possibilities for a node: first, with probability 



t+1 



it may get an extra edge from the new node while its ex- 
isting edges remain undeleted, and thus increase its own 
degree by 1; and second, with the complimentary prob- 
ability 1 — tiri , the nodes may remain in the former 
state with the former degree. It should be noted that Eq. 
(fl]) and all the following ones are exact for alH > 2. 

The total degree distribution of the entire network can 
be obtained as 



P(k,t) = T ^Y /P (k,s,t) 



(2) 



Using this and applying X)*=o ^° both sides of Eq. fl} , 
we get the following master equation for the degree dis- 
tribution: 

(t + 2)P(k,t + l)-{t + l)P(k,t) 
= 2(1 - q)P(k - 1, t) - 2(1 - q)P(k, t) + 4,2- (3) 

The corresponding stationary equation, i.e., at t — > oo, 
takes the form 



III. TOPOLOGICAL PROPERTIES 

We focus on the behavior of the topological character- 
istics, in terms of the degree distribution, the clustering 
coefficient, and the average path length, as a function of 
the parameter q. 



A. Degree Distribution 

The degree distribution is one of the most important 
statistical characteristics of a network. For q = 1, all 
nodes have the same number of connections 2, the net- 
work exhibits a completely homogeneous degree distri- 
bution. Next we focus the case Q < q < 1. In or- 
der to conveniently describe the computation of the net- 
work characteristics, we label nodes by their birth times, 
s = 0, 1, 2, . . . , t, and use p(k, s, t) to denote the probabil- 
ity that at time t a node created at time s has a degree 
k. At time t, there are t+1 internode intervals along 
the ring and each node has two intervals (one on either 
side) . The master equation (25l . [26| governing the evolu- 
tion of the degree distribution of an individual node has 



(3 - 2q)P(k) - (2 - 2q)P(k - 1) = 5 ka . (4) 

Eq. ((3]) implies that P(k) is the solution of the recursive 
equation 

P(k) = ll^Yq^ k -^ fOTfc>2 (5) 
[ 1/(3 -2q) forfc = 2 

giving 

which decays exponentially with k. For q = 0, Eq. (p]) 
recovers the result previously obtained in [l4]]. Thus the 
resulting network is an exponential network. Note that 
most small-world networks including the WS model be- 
long to this class [1 [3 [H [3 . 

In Fig. [21 we report the simulation results of the degree 
distribution for several values of q. From Fig. [21 we can 
see that the degree spectrum of the networks is continu- 
ous and the degree distribution decays exponentially for 
large degree values, in agreement with the analytical re- 
sults and supporting a relatively homogeneous topology 
similar to most small- world networks [TJ, [IH, [l|| ■ 
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FIG. 2: Semilogarithmic graph of degree distribution of the 
networks with order N = 10 5 . The solid lines are the analytic 
calculation values given by Eq. (|6}. 



B. Clustering coefficient 

Most real-life networks show a cluster structure which 
can be quantified by the clustering coefficient [l|, H, [3, 0, 
5] . The clustering coefficient of a node gives the relation 
of connections of the neighborhood nodes connected to 
it. By definition, clustering coefficient Ci of a node i is 
the ratio of the total number a of existing edges between 
all ki its nearest neighbors and the number ki [ki — 1)/2 of 
all possible edges between them, ie Ci = 2ei/[ki(ki — 1)]. 
The clustering coefficient C of the whole network is the 
average of all individual C^s. 

For the case of q — 1, the network is a one-dimensional 
chain, the clustering coefficient of an arbitrary node and 
their average value are both zero. 

For the case of q = 0, using the connection rules, it is 
straightforward to calculate exactly the clustering coeffi- 
cient of an arbitrary node and the average value for the 
network. When a node i enters the network, ki and 
are 2 and 1, respectively. After that, if the degree ki in- 
creases by one, then its new neighbor must connect one of 
its existing neighbors, i.e. increases by one at the same 
time. Therefore, is equal to ki — 1 for all nodes at all 
time steps. So there exists a one-to-one correspondence 
between the degree of a node and its clustering. For a 
node v with degree k, the exact expression for its cluster- 
ing coefficient is -| , which has been also been obtained in 



other models [hUE [M [27|, [l| . This expression for the 
local clustering shows the same inverse proportionality 
with the degree as those observed in a variety of real- 
life networks |29j. In this limiting case, the clustering 
coefficient C of the whole network is given by 



oo 



0.6479. 



(7) 
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So in the limit of large t the clustering coefficient is very 



FIG. 3: The clustering coefficient C of the whole network as 
a function of q. 



high. 

In the range < q < 1, it is difficult to derive an ana- 
lytical expression for the clustering coefficient either for 
an arbitrary node or for the average of them. In order 
to obtain the result of the clustering coefficient C of the 
whole network, we have performed extensive numerical 
simulations for the full range of q between and 1. Sim- 
ulations were performed for system sizes 10 5 , averaging 
over 20 network samples for each value of q. 

In Fig. [3l we plot the clustering coefficient C as a func- 
tion of q. It is obvious that C decreases continuously 
with increasing q. As q increases from to 1, C drops 
almost linearly from 0.6479 to 0. Note that although the 
clustering coefficient C changes linearly for all q, we will 
show below that in the large limit of q, the average path 
length changes exponentially as q . This is little different 
from the phenomenon observed in the WS model where 
C remains practically unchanged in the process of the 
network transition to a small world. 



C. Average Path Length 

Certainly, the most important property of a small- 
world network is a logarithmic average path length 
(APL) (with the number of nodes). Here APL means the 
minimum number of edges connecting a pair of nodes, av- 
eraged over all pairs of nodes. It has obvious implications 
for the dynamics of processes taking place on networks. 
Therefore, its study has attracted much attention. 

For the case of q = 1, the average path length increases 
linearly with network size. For the case of q = 0, the 
network grows stochastically. Generally speaking, for a 
randomly growing network, the analytical calculation for 
APL is difficult. Below, we will give an upper bound 
for the APL of this particular case, which shows that the 
APL increases at most logarithmically with network size. 

If d(i,j) denotes the distance between nodes i and j, 
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we introduce the total distance of the network with size 
N as a(N): 



3000 



a(N)= 

0<i<j<N-l 

and we denote the APL by L(N), defined as: 

2a(N) 



L(N) = 



N(N - 1)' 



(8) 



(9) 



In the limiting case q — 0, the distances between existing 
node pairs will not be affected by the addition of new 
nodes. Then we have the following equation: 



N-l 



a(N + l) = a(N)+J2 d(i,N). 



(10) 



Assume that the node N is added and connected to two 
nodes w\,W2 linked by edge E, then Eq. ([10]) can be 
rewritten as: 



JV-l 



a(N + l) = cr(N)+J2{D(i,w) + l] 

i=0 

N-l 

= a(N)+N+ Y,D(i,w), (11) 



i=0 



where D(i, w) = min{d(i, tui), <2(i, i/^)}- Constricting the 
edge E continuously into a single vertex w (here we as- 
sume that w = wi), we have D(i,w) = d(i,w). Since 
d(wi, w) = d{w2, w) = 0, Eq. pip can be rewritten as: 



a(N + 1) = a(N) + N + ^ d(i, w), 



(12) 



where T = {0, 1, 2, • • • , N — 1} — {wi, w 2 } is a node set 
with cardinality N — 2. The sum d(i,w) can be 

considered as the total distance from one node w to all 
the other nodes in the network with size N — 1, which 
can be roughly evaluated by mean-field approximation in 
terms of L(N - 1) as 0, IM HI : 



Y,d(i,w)n(N-2)L(N-l). 



(13) 



ter 



Note that, as L(N) increases monotonously with N, it is 
clear that: 



x / x 2a(N-l) 2a(N) 
(N — 2)L(N — 1) = — ^ — ^< v ; 



N-l 



N 



(14) 



Combining Eqs. (fT2|) . (fT3|) and (fT4|) . one can obtain the 
inequation: 



a(N + 1) < a(N) +N + 



2a(N) 
N ' 



(15) 



Considering Eq. (|15p as an equation and not an inequal- 
ity, we can provide an upper for the variation of cr(N) 
as 



da(N) 
dN 



N + 



2a(N) 



(16) 
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FIG. 4: Average path length L versus parameter q. 



which leads to 



a(N) = N 2 (lnN + a), 



(17) 



where a is a constant. As a(N) ~ A^ 2 ^^, we have 
L(N) - In N. Note that as we have deduced Eq. 1(17) 
from an inequality, then L(N) increases at most as IniV 
with N. Therefore, we have proved that in the special 
case of q = 0, there is a slow growth of APL with network 
size N. 

For < q < 1, in order to obtain the variation of 
the average path length with the parameter q, we have 
performed extensive numerical simulations for different q 
between and 1. Simulations were performed for system 
sizes 10 4 , averaging over 20 network samples for each 
value of q. In Fig. [H we plot the average path length 
L as a function of q. We observe that, when lessening 
q from 1 to 0, average path length L drops drastically 
from a very high value a small one, which predicts that a 
phase transition from large- world to small- world occurs. 
This behavior is similar to that in the WS model. 

Why is the average path length L low for small ql The 
explanation is as follows. The older nodes that had once 
been nearest neighbors along the ring are pushed apart as 
new nodes are positioned in the interval between them. 
From Fig. [1] we can see that when new nodes enter into 
the networks, the original nodes are not near but, rather, 
have many newer nodes inserted between them. When q 
is small, the network growth creates enough "shortcuts" 
(i.e. long-range edges) attached to old nodes, which join 
remote nodes along the ring one another as in the WS 
model 0]. These shortcuts drastically reduces the aver- 
age path length, leading to a small-world behavior. 



IV. CONCLUSIONS 

In summary, we have proposed a one-parameter model 
of growing small-world networks. In our model, in ad- 
dition to new edges connecting new nodes and old ones, 
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edges between old nodes may be removed. The presented 
model interpolates between one-dimensional regular ring 
and small-world networks, which allow us to explore the 
crossover between the two limiting cases. We have ob- 
tained both analytically and numerically the solution for 
relevant parameters of the network and observed that our 
model exhibits the classical phenomenon as that in the 
WS model. Our model may provide a useful tool to inves- 
tigate the influence of the clustering coefficient or average 
path length in different dynamics processes taking place 
on networks. In addition, using the idea presented here, 



one can also construct models interpolating between ho- 
mogeneous and heterogeneous networks 33] . 
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